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Abstract 

o 

Experimental results for sub-barrier nuclear fusion reactions show cross sec- 
tion enhancements with respect to bare nuclei which are generally larger than 
those expected according to electron screening calculations. We point out that 
energy spread of target or projectile nuclei is a mechanism which generally 
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CO . provides fusion enhancement. We present a general formula for calculating 
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the enhancement factor and we provide quantitative estimate for effects due 
to thermal motion, vibrations inside atomic, molecular or crystal system, and 
due to finite beam energy width. All these effects are marginal at the energies 
which are presently measurable, however they have to be considered in future 
experiments at still lower energies. This study allows to exclude several ef- 
fects as possible explanation of the observed anomalous fusion enhancements, 
?-) | which remain a mistery. 

I. INTRODUCTION 

x. 

The chemical elements were created by nuclear fusion reactions in the hot interiors of 
remote and long- vanished stars over many billions of years [1]. Thus, nuclear reaction rates 
are at the heart of nuclear astrophysics: they influence sensitively the nucleosynthesis of the 
elements in the earliest stages of the universe and in all the objects formed thereafter, and 
they control the associated energy generation, neutrino luminosity and evolution of stars. A 
good knowledge of their rates is essential for understanding this broad picture. 

Nuclear reactions in static stellar burning phases occur at energies far below the Coulomb 
barrier. Due to the steep drop of the cross section a(E) at sub-barrier energies, it becomes 
increasingly difficult to measure it as the energy E is lowered. Generally, stellar fusion rates 
are obtained by extrapolating laboratory data taken at energies significantly larger than 
those relevant to stellar interiors. Obviously such an "extrapolation into the unknown" 
can lead to considerable uncertainty. In the last twenty years a significant effort has been 
devoted to the experimental exploration of the lowest energies and new approaches have 
been developed so as to reduce the uncertainties in the extrapolations. In particular, the 
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installation of an accelerator facility in the underground laboratory at LNGS [2] has allowed 
the cr(E) measurement of 3 He( 3 He, 2p) 4 He down to its solar Gamow peak, E ±A/2 = (21±5) 
keV [3] so that for this reaction no extrapolation is needed anymore. 

As experiments have moved well down into the sub-barrier region, the screening effect 
of atomic electrons has become relevant [4-6]. With respect to the bare nuclei case, the 
Coulomb repulsion is diminished, the tunneling distance R t is reduced, and the fusion prob- 
ability, which depends exponentially on R t , is enhanced. The electron effect on the reaction 
can be seen as a transfer of energy U (the screening potential energy) from the electronic 
to the translational degrees of freedom. For each collision energy E, one has an effective 
energy E eB = E + U and a cross section enhancement: 

f = a(E + U)/a(E) (1) 

The screening potential energy U is easily estimated in two limiting cases [5]: In the 
sudden limit, when the relative velocity v re i of the nuclei is larger than the typical electron 
velocity t>o = e 2 /h: the electron wave function during the nuclear collision is frozen at the 
initial value \l/i n and the energy transferred from electrons to the nuclei is thus 

U m = {* in \Z 1 e 2 /r u \*in) , (2) 

where here and in the following the index 1 (2) denotes the projectile (target) nucleus and a 
sum over the electrons is understood. In the adiabatic limit, i.e. when v re i <C t>o: electrons 
follow adiabatically the nuclear motion and at any internuclear distance the electron wave 
function \l/ a( j corresponds to an energy eigenstate calculated for fixed nuclei. As the nuclei 
approach distances smaller than each atomic radius , \l/ a d tends to the united atom (i.e. with 
nuclear charge Z = Z x + Z 2 ) limit, $ un . The kinetic energy gained by the colliding nuclei is 
thus 

U^d = ^in e un i (3) 

where e in (e un ) is the electron energy of the isolated (united) atom in the corresponding 
states. 

We like to stress a few important features: 

i) Screening potential energies, which are in the range 10-100 eV, are definitely smaller 
than the practical collision energies (1-100 keV), nevertheless they can produce appreciable 
fusion enhancements due to the exponential dependence of the cross section. 

ii) In the adiabatic limit the electron energy assumes the lowest value consistent with 
quantum mechanics. Due to energy conservation, the energy transfer to the nuclear motion 
is thus maximal in this case (U < U a a) and the observed cross section enhancement should 
not exceed that calculated by using the adiabatic potential: 

f<f a d = <r(E + U ad )/a(E) (4) 

in) The enhancement factors which have been measured are generally larger than ex- 
pected. A summary of the available results is presented in Table I. The general trend is 
that the enhancement factors exceed the adiabatic limit. Recent measurement of d(d, p) 3 H 
with deuterium implanted in metals [7] have shown enhancements of the cross sections with 
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respect to the bare nuclei case by factors of order unity, whereas one expects a few percent 
effect. In other words, if one derives an "experimental" potential energy U cx from a fit of ex- 
perimental data according to eq. (1), the resulting values significantly exceed the adiabatic 
limit U ad . In the case of deuterium implanted in metals, values as high as U ex ~ 700 eV have 
been found [7], at least an order of magnitude larger than the expected atomic value U a( i- 
Several theoretical investigations have resulted in a better understanding of small effects in 
low energy nuclear reactions, but have not provided an explanation of this puzzling picture. 

iv) Dynamical calculations of electron screening for finite values of the relative velocity 
show a smooth interpolation between the extreme adiabatic and sudden limits [8,9]. In 
fact, one cannot exceed the value obtained in the adiabatic approximation because the 
dynamical calculation includes atomic excitations which reduce the energy transferred from 
the electronic binding to the relative motion. 

v) The effects of vacuum polarization [10], relativity, Bremmstrahlung and atomic po- 
larization [11] have been studied. Vacuum polarization becomes relevant when the minimal 
approach distance is close to the electron Compton wavelength but it has an anti-screening 
effect, corresponding to the fact that in QED the effective charge increases at short distances. 
All these effects cannot account for the anomalous enhancements. 

Although one cannot exclude some experimental effect, e.g. a (systematic) overestimate 
of the stopping power, the general trend is that most reactions exhibit an anomalous high 
enhanchement. Phenomeno logically, this corresponds to an unexplained collision energy 
increase in the range of 100 eV. 

Actually, the anomalous experimental values U ex look too large to be related with atomic, 
molecular or crystal energies. Some other processes, involving the much smaller energies 
available in the target, should mimic the large experimental values of U. As an example, if 
the projectile approaches a target nucleus which is moving against it with energy E 2 <C E, 
the collision energy is increased by an amount: 



For d + d reactions at (nominal) collision energy E — 10 keV, a target energy E 2 = 0.5 eV 
is sufficient for producing U = 100 eV. 

Generally, one expects that opposite motions of the target nuclei are equally possible. 
Even in this case, however, the effect is not washed out: due to the strong nonlinearity of 
the fusion cross section the reaction probability is much larger for those nuclei which are 
moving against the projectile. 

In this spirit, we shall consider processes associated with the energy spread of the colliding 
nuclei. These processes generally lead to an enhancement of the fusion rate, for the reasons 
just outlined. 

In the next section we shall first consider the thermal motion of the target nuclei. For this 
example, we shall derive an expression for the enhancement factor on physical grounds and 
then we shall outline the effects of an energy spread for the extraction of the astrophysical 
S-factor from experimental data. 

The treatment is generalized in sect. Ill and in sect. IV it is applied to study energy 
spreads due to motion of the nuclei inside atoms, molecules and crystals. Beam energy width 
and straggling are also considered. 




(5) 
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In summary, all the effects turn out to be too tiny to explain the observed anomalous 
enhancements. Nevertheless, they have to be considered in analyzing the data, particularly 
in future experiments at still lower energies. 



II. THE EFFECT OF THERMAL MOTION OF TARGET NUCLEI 

In this section we consider the effects of thermal motion of the target nuclei. We shall 
make several simplifications, in order to elucidate the main physical ingredients. In this way 
we shall derive a simple expression for the enhancement factor on physical grounds. 

Essentially, we shall concentrate on the exponential factor of the fusion cross section, 
neglecting the energy dependence of the pre-exponential factors, and we shall only consider 
the effect of the target motion in the direction of the incoming particle, neglecting the trans- 
verse motion. When these simplifications are removed the result is essentially confirmed: 
see the more general treatment of sect. III. 

The fusion cross section at energies well below the Coulomb barrier is generally written 

as: 



where E = \^v 2 el is the collision energy, ji = mim 2 /(mi + m 2 ) is the reduced mass, Vo = 
ZiZ 2 e 2 /h and S(E) is the astrophysical S-factor 1 . The cross section is more conveniently 
expressed in terms of the relative velocity of the colliding nuclei i> re i: 



At energies well below the Coulomb barrier, v re \ <C Vq, the main dependence is through the 
exponential factor, so we shall treat the pre-exponential term as a constant: 



We consider a projectile nucleus with fixed velocity V impinging against a target where 
the nuclei have a thermal distribution of velocity. Since the target nucleus velocity v is 
generally much smaller than V = |V|, one can expand l/i> re i — 1 / — v| and retain the 
first non vanishing term: 




(6) 




(7) 




(8) 




(9) 



where v\\ is the target velocity projection over the V-direction. 



For convenience of the reader, we recall that vq = e 2 /fi and thus Vq = Z1Z2VQ. 
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The enhancement factor with respect to the fixed target case, / = (a)/a(V), is thus 
calculated by averaging exp(— V v \\/V 2 ) over the v\\ distribution: 



P("ll)= ,, expl-i^l, (10) 



where (v 2 ) = kT/m 2 . The integral 



CO 




is easily evaluated by using a (saddle point) trick similar to that used by Gamow for evaluat- 
ing stellar burning rates. The product of the Gaussian and the exponential functions (Fig.l) 
results in an (approximately) Gaussian with the same width, centered at vq = — (vpVo/V 2 , 
its height giving the enhancement factor: 

f = exp [-w^) ■ (12) 

Concerning this equation, which is the main result of the paper, several comments are 
needed: 

%) Since the term in parenthesis in eq. (12) is positive, one has / > 1, i.e. the energy 
spread always results in a cross section enhancement. One cannot ignore the target velocity 
distribution for the calculation of the reaction yield since nuclei moving towards the projectile 
have a larger weight in the cross section. 

ii) The main contribution to the cross section comes from target nuclei with velocity 
close to vq. When V < \J {v 2 ) 1 / 2 Vo, this velocity is larger than the typical thermal velocity 
(vjj) 1//2 . This result is equivalent to the Gamow peak energy in stars, which is significantly 
higher than the thermal energy kT. In terms of the energy, by putting E 2 = \raiVQ in 
eq. (5), we see that the "most probable" collision energy is 2 : 

E„^E + 2(^-)^E t , (13) 
Vmi + m 2 J V 

where E t = |m 2 (f 2 ) = \kT is the average thermal energy associated with the motion in the 
collisional direction. 



in) The energy dependence of eq. (12), 



/ = exp 





( ml ) 


E t E ~ 




Vmi + m 2 J 


E 2 \ 



(14) 



where Eq = l/^V 2 , is different from that resulting from electron screening / = exp(D/E 3 ^ 2 ). 



2 The most probable energy E mp has not to be confused with the effective energy E e ff. 
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iv) The resulting effects are anyhow extremely tiny. For example, for d + d collisions 
(V = e 2 /h) at E — 1 keV (V = 1/5 Vo) and room temperature ((v 2 ) 1 ^ 2 = 5 ■ 10~ 4 V ) one 
has / — 1 ~ 10~ 4 . A 10% enhancement would correspond to kT ~ 30 eV. 

v) The same method can be extended to other motions of the target nuclei, provided 
that the velocity distribution is approximately Gaussian and if other interactions of the 
nuclei during the collision are neglected (sudden approximation). One has to replace (v?<) in 
eq. (12) with the appropriate average velocity associated with the motion under investigation. 
Vibrations of the target nucleus inside a molecule or a crystal lattice can be treated in this 
way, since the vibrational times are much longer than the collision times. These and other 
similar effects will be discussed in sect. IV. 

vi) From the discussion presented above one gets an easy procedure to correct the 
experimental results for taking into account the effect of an energy spread. If the as- 
trophysical S-factor has been measured at a nominal collision energy E = ^pV 2 , from 
S'cxp = cr cxp E exp(Vo/V), then the "true" S-factor is obtained as S = S cxp /f, where / is 
given by eq. (12) and the "true" energy is changed from E to E mp given in eq. (13) ( Fig. 2). 
In summary, the effect of the energy spread translates into both a cross section enhancement 
and an energy enhancement. 



In this section we shall provide a more general discussion of the energy spread effects, 
which will substantially confirm eq. (12) and which can be applied to a rather large class of 
processes. The main assumption is that the projectile motion is fast in comparison with the 
other motions, so that the sudden approximation can be used. 

Let us consider a projectile with velocity V impinging onto a thin target (density n 
and thickness L), where energy loss can be neglected. The interaction probability P is the 
product of the interaction probability per unit time p = n(av re \) with the time spent in the 
target, L/V. The measured counting rate A = elp, where / is the beam current and e is the 
detector efficiency, is thus: 



As in stars, the quantity which is physically relevant is thus (av TC \), where the average has 
to be taken over the target nuclei velocity distribution. 

This distribution is due to the coupling with other degrees of freedom. Inside an atom 
(or a molecule, or a crystal) the nucleus is vibrating, its motion is altered by the arrival of 
the projectile nucleus and the calculation of the average is complicated in the general case. 
However, if the velocity V of the impinging particle is large in comparison with the velocity 
v of the target nucleus, the problem is simplified. The target wave function does not have 
time for significant evolution during the collision and it can be taken as that of the initial 
(unperturbed) state. This is the main content of the sudden approximation: the velocity 
distribution of the target nuclei p{v) can be taken as the initial one pi n (v) and one has to 
compute: 



III. GENERAL TREATMENT 




(15) 
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A = 



IenL 
V 



d 3 v p in (v) a(v Te i)v 



rel 



By using eq.(7), one has thus to compute: 

IenL 



A = J d 3 v p in (v) 



2S(v rcl ) ( V 
exp 

pV Te \ \ V re \ 



(16) 



(17) 



We recall that S is a weakly varying function of energy so that it can be taken out of the 
integral. 

Since we are assuming V 2 ^> (v 2 ), we expand the integrand g — ^ exp(— Vo/i> re i) in 
powers of v and keep the lowest order terms: 



A 



2SIenL 
pV 



J d 3 v p in (v) 



g v =o + Vi(dig) v=0 + -ViVj{didjg) v=0 



(18) 



We shall consider distributions which are symmetrical for inversions and rotations around 
the collision axis V. In this case the term linear in v vanishes and the result is: 



A = 



ISlenL 



pV 2 



, v 

exp ( -- 



(19) 



{1 + 



<^o 2 



2 V 4 



(vl)Vo 2 
2V A 



V_ 

V 



4— + 2 — 

V Wo 

Vo 



} , 



where the index || (±) denotes the component of the velocity along (transverse to) the 
collision axis. 

The term in front of the curl bracket is the counting rate calculated neglecting the target 
energy spread. So, if we define the enhancement factor / as the ratio of the measured 
counting rate A to the rate calculated for fixed velocity Ay, we have: 



d 3 v p in (v) 



1 

v re \ 



exp 



we have now: 



2 V 4 

(vl) Vq 2 
2 V 4 



V 

1 -4— + 2 
Vo 



V_ 

V 



Vo 



+ 



(20) 



(21) 



V_ 

Vo 



For a one dimensional motion (v± = 0) it simplifies to: 



7 = 1 + 



2V 4 



2 r 



V 



21 



(22) 



For the case of a spherically symmetrical distribution, (v 2 ) = l/2{v 2 _) one gets: 
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/ = 1 + 



(vj) 
2 V 4 



(23) 



This equation can be easily compared with the result of the previous section concerning the 
thermal energy effect. By expanding eq. (12) one gets: 



,_ 1 , l(j)^o 2 
/ 2 V 4 



(24) 



This is the same as eq. (23) apart for the last term which is negligible at small velocities, 
since it is a higher order contribution in V/Vq. Note that this last term arises from the 
variation of the pre-exponential factor l/i> re i, which was neglected in the simplified treatment 
of sect. II. Clearly this term, once averaged over the target distribution, is smaller than 1/V 
and therefore it provides a reduction of the rate, as implied by the negative coefficient in 
eq.(23). 

The previous results have been obtained by neglecting higher order terms in the expansion 
of g. Their contribution is suppressed by a factor (vf) V 2 /V 4 .Thus the previous results are 

not valid for V <$C 



(v 2 ) 1 / 2 V , as can be simply understood. In this case, one cannot expand 
the integrand function g(v), since it changes faster than the distribution function p(v) over 
a large range of target velocities. More precisely, the decrease of p(v) is counter-balanced 
by the increase of g(v) in a velocity range which is typically larger than the average target 
velocity dispersion (v 2 ) 1 ^ 2 . As a consequence, the tails of the distribution function p{y) give 
a relevant contribution to the counting rate, leading to an increase of the factor / with 
respect to the simple estimate eq.(23). 

It is difficult to obtain a general expression for / in this low velocity regime. The factor 
/ depends, in fact, on the shape of the distribution function. In the case of a gaussian 
distribution function, p(v) oc exp(— v 2 /2(v 2 )), one can use the Gamow "trick" described in 
the previous section which leads to eq.(12). For distribution functions which decrease more 
slowly with v one expects larger effects. 

In order to have, however, a general result for the low velocity ((v 2 ) < V 2 < {v 2 ) 1 ^ 2 Vq) 
behaviour of /, we note that, being the counting rate A an increasing function of the pro- 
jectile velocity V, one has: 



A > A(0) 



IenL 
V 



J d 3 v pi, 



1 



cxp 



V 



This means that the enhancement factor / should be larger than: 

■Vn 



f = Vexp (^yj J d 3 v p in {v 



1 / V 
- exp 



v 



V 



(25) 



(26) 



IV. APPLICATIONS 

The method developed in the previous sections, summarized in eq. (12) or in the more 
accurate eq. (21), can be applied to several motions of the target nuclei (vibrations inside 
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an atomic, molecular or crystal system), provided that interactions with other degrees of 
freedom during the collision can be neglected. Simply, one has to compute the value of (v 2 ) 
which is appropriate to the system under consideration. Also, the treatment can be easily 
extended to the effect of beam energy width and straggling. 



P-M = -S ,..2 ,°,2u > ( 28 ) 



A. Nuclear motion inside the atom 

Very much as the motion of a star in the sky is affected by the presence of planets 
around it, the nucleus inside an atom is vibrating around the center of mass of the atomic 
system. The nuclear momentum distribution P(p) is immediately determined from that of 
the atomic electrons P e (p e ) by requiring that the total momentum of the atom vanishes in the 
center of mass (p = —p e ), where p e is the (total) momentum carried by the electron(s), i.e. 
P{p) = P e (—p e ) and the initial nuclear velocity distribution p m {v) is immediately determined 
from v = p/m.2, where m 2 is the target nucleus mass. 

For the case of Hydrogen (isotope) in the ground state, the atomic electron momentum 
distribution is: 

P(v) = - {meVo)5 (27) 

so that the nucleus velocity distribution is: 

_8_ 

7T 2 {v 2 +u 2 y 

where uq = (m e /m 2 )t>o = (m e /m 2 )e 2 //i, is the typical velocity associated with the target 
nuclear motion. In practice, this is definitely smaller than the collision velocity V, so that 
the sudden approximation holds and the results of the previous section can be applied. 
One can easily evaluate that: 

«>4-H© 2 "» ■ (29) 

so that for Hydrogen- Hydrogen (or deuterium-deuterium) collisions, for which Vq = 
Z 1 Z 2 e 2 /h = Vq, by using eq. (23), one obtains for the enhancement factor : 

6 \m 2 ) \V ) \ V 

This is an extremely tiny correction, since one has / at — 1 — 2 • 10~ 5 for a d-d collision at 
E — 1 keV energy . 

In the low energy regime, i.e. when V < V M o Vo — (m e /m 2 ) 1 ^ 2 v , the previous estimate 
has to be corrected to take into account the contribution of the tails of the distribution 
function. By using eq.(26) we can easily estimate: 

In Fig. 3 we compare the approximate expressions with the numerical evaluation of eq. (21). 
In the whole range a good approximation to the full numerical calculation is provided by 
/ — fat + (/at)o- 



f« = l + U~) (tf) f 1 " 2 ^) • (30) 
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B. Molecular vibrations 



Let us consider, as an example, reactions involving a deuterium nucleus bound in a D 2 
molecule. The target nucleus is vibrating, the vibration energy in the ground state being 
-Evib = 0.19 eV . This energy is shared between the two nuclei and between potential and 

l/4£ vib . 



kinetic energy, so that the average kinetic energy of each nucleus is |m d (t> 2 ) vib 



The target nucleus velocity, (v ) 



vib 



10 



-6 2 
y 0> 



is much smaller than the projectile velocity 



so that the sudden approximation applies again. By using eq. (12) and assuming a random 
orientation of the molecular axis , (v?<) = l/3(t> 2 ) v ib, we get: 



mol 



cxp 



QV 4 



(32) 



This corresponds to a 10 4 correction at E=l keV. Conversely, an enhancement correction 
of 10% would correspond to E vib ~ 200 eV. 



C. Local vibrations in a crystal lattice 

When a deuterium nucleus is implanted in a crystal, it generally occupies an interstitial 
site where it performs local vibrations. The vibration energy E cr depends on the host lattice, 
being typically in the range of 0.1 eV, very similar to the molecular vibration scale. Effects 
associated with vibrations in the crystal are thus similar to those calculated for the D 2 
molecule: 

fcr — /mol • (33) 



D. Finite beam width and straggling 



In an ideal accelerator all projectiles have the same energy E^b- Actually, due to several 
physical processes (voltage fluctuations, different orbits...) the beam will have a finite energy 
width A. As an example, in the LUNA accelerator one has A ~ 10 eV. Furthermore, 
when the beam passes through the target, fluctuations in the energy loss will produce an 
enlargement of the energy width (straggling). Thus, even neglecting the target motion, there 
is a collision energy spread. The beam energy distribution, 



P(E') ~ exp 



(E — EiabY 

2A 2 



(34) 



gives a velocity ditribution with: 



m d E } 



lab 



(35) 
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By using eq. (21) the enhancement factor is thus 3 : 



/ = exp 



U 2 A 2 



(36) 



Effects are very small in the case of LUNA: for d+d at E — 1 keV and A = 10 eV one 
has / — 1 ~ 2 • 1CT 5 . The effect behaves quadratically with A and it can be significant if 
momentum resolution is worse. Conversely, an enhancement correction of 10% corresponds 
to A ~ 250 eV. 



E. Polynomial velocity distributions 

One could suspect that velocity distributions of different shape can provide enhancements 
significantly larger than the tiny effects which we have found so far. 

In this spirit, let us consider the case of a polynomial velocity distribution, 

M = {V 2 f Br (37) 

where the slowly decreasing tail should provide a significant enhancement. Clearly the more 
favourable cases correspond to small values of n. The requirement that (v 2 ) is finite implies 
n > 3, so we consider n = 3 in order to maximize the tail effect. The normalized distribution 
is in this case: 



P{v) = ^¥i 2 [v 2 + {imv 2 )f (38) 



The low energy enhancement factor / of eq. (26) becomes now: 

In order to have / ~ 1.1 for d+d collisions at E — 1 keV one needs (v 2 ) ~ 3 • 10~ 2 V 2 , 
which corresponds to an average energy in the range of 1 keV, well above the physical scale 
of the process. 



V. CONCLUDING REMARKS 

We summarize the main points of this paper: 



3 For the sake of precision, the counting rate is now: A = eInL(a)b earn . This is different from eq. 
(15). A calculation of the average, similar to that presented in sect. Ill, yields the same expression 
as in eq.(19) for the leading term in V/Vq and different numerical coefficients for the higher order 
(negligible) terms. 
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i) Energy spread is a mechanism which generally provides fusion enhancement. 
ii) We have found a general expression for calculating the enhancement factor /: 



/ = exp 



' Z 1 Z 2 e 2 \ 2 (v l{ ) 



h 



2V A 



(40) 



Hi) We have provided quantitative estimates for the enhancement effects. For a d+d 
collision one has: 

thermal motion: / — 1 ~ 1CT 4 (E/lkeV)~ 2 
vibrational motion: / - 1 ~ (1(T 5 - 1(T 4 ) (E/lkeV)" 2 
beam width: / - 1 ~ (E/lkeV)" 3 

iv) All these effects are marginal at the energies which are presently measurable, however 
they have to be considered in future experiments at still lower energies. 

v) This study allows to exclude several effects as possible explanation of the observed 
anomalous fusion enhancements, which remain a mistery. 
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TABLES 



TABLE I. Summary of results for electron screening effects 



Reaction 


Uex (eV) 


U ad a (eV) 


Ref. 


d(d,p)t 


25 ±5 b 


28.5 


[12] 


3 He(d,p) 4 He 


219 ±7 


114 


[13] 


d( 3 He,p) 4 He 


109 ±9 


102 


[13] 


3 He( 3 He,2p) 4 He 


294 ± 47 


240 


[3] 


3 He( 3 He,2p) 4 He 


432 ± 29 


240 


[14] 


6 Li(p,a) 3 He 


470 ± 150 


184 


[15] 


6 Li(d,a) 4 He 


320 ± 50 


184 


[16] 


7 Li(p,a) 4 He 


330 ± 40 


184 


[17] 


9 Be(p,d) 8 Be 


900 ± 50 


262 


[18] 


11 B(p,a) 8 Be 


430 ± 90 


346 


[19] 


a Values calculated for atomic target, following [5]. It is assumed that at fusion hydrogen projectiles are 
charged or neutral with equal probability. Helium projectiles are assumed to be He + (He) with 20% (80%) 
probability. 

6 This value results from gaseous target. Much larger values have been found when deuterium is implanted 



in metals [7]. 
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FIGURES 




v K /<v K 2 > 1 



FIG. 1. A sketch of the contribution to the averaged cross section. p(v\\) is defined in eq.(10) 
and Q = exp(—vo(v\\}/V 2 ). 
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FIG. 2. Extraction of the S*-factor from experimental data. 
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FIG. 3. Fusion enhancement due to nuclear motion inside a H atom, we present the numerical 
evaluation of eq.(21) (full line), the approximations of eq.(12) (dotdashed) and of eq.(23) (dotted), 
the low velocity limit eq.(31) (dashed). 
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